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Abstract. We prove existence, uniqueness and asymptotics of global smooth 
solutions for the Landau-Lifshitz-Gilbert equation in dimension n > 3, valid 
under a smallness condition of initial gradients in the L" norm. The argu- 
ment is based on the method of moving frames that produces a covariant 
complex Ginzburg-Landau equation, and a priori estimates that we obtain by 
the method of weighted-in-time norms as introduced by Fujita and Kato. 

Introduction and results 

The Landau-Lifshitz-Gilbert equation (LLG) is the fundamental evolution law in 
magnetism. It governs the dynamics of continuous spin systems, i.e. of director 
fields m = (0, oo) x M" — >• with values in the unit sphere §^ C M'^ where 
typically n — 3. In its original form as introduces by Landau and Lifshitz [231 113] 
the equation reads 

dm , ^ ^ 

(1) — — = m X Am — Am X m X Am 
^ ' dt 

where A > is a damping parameter and x is the vector product on M.^ . Observe 
that — m X m X Am = Am -I- |Vmpm, which is sometimes called the tension field. 
Mathematically LLG can be considered as a dissipative version of the Schrodinger 
flow for harmonic maps into (where A = 0). Emphasizing its parabolic char- 
acter, LLG can also be considered as a quasilinear perturbation of the heat flow 
for harmonic maps by the (conservative) precession term — m x Am. A common 
Liapunov functional is given by the Dirichlet energy 



E{m) = - / iVmpda;. 



In dimension n = 2, LLG is energy critical, which means that the scaling symmetry 
m,^{t,x) = m{t/e'^,x/e) for e > featured by ([T|) preserves the Dirichlet energy, 
i.e. E{mg{t)) = E{m{t/ e^)) for all t > {). Global existence of weak solutions and 
uniqueness in the class of energy decreasing solutions can be shown with arguments 
parallel to those for the heat flow, see [TH [TJl [TB] . In contrast, the existence of finite 
time singularities as known for the heat flow has not been shown conclusively, but 
there is at least strong numerical evidence, see [2]. The possible blow-up cenario, 
however, is precisely characterized through the development of bubbles at energy 
concentration points. Consequently, initial energies below 47r, which is the topolog- 
ical lower bound for a full cover of S^, will guarantee global regularity. Large energy 
solutions on the other hand are particularly interesting if topologically nontrivial 
data, featuring e.g. magnetic vortices, is concerned. Recent results show that in the 
presence of strong potentials of Ginzburg-Landau type, conventional energy bounds 
can be bypassed by certain well-preparedness assumptions on the initial data that 
prevent the formation of (extra) bubbles, see [1TJ[5D]. 
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In dimension n > 3, LLG becomes super-critical with respect to the Dirichlet en- 
ergy, i.e. with as before E{ms{t)) — J' as e — ?> 0. Whereas the proof of existence 
of weak solutions remains unaffected, see [1] , questions of uniqueness and regularity 
become more delicate. In particular considerations solely based on the energy are 
insufficient for ruling out concentration effects, and new adapted quantities have to 
come into play. More advanced features of the harmonic map heat flow equation, 
such as Bochner identities and monotonicity formulas giving rise to crucial a priori 
estimates, see [3T], are not at hand for LLG. In dimension n = 3,4 perturbative 
methods based on elliptic problems on suitable time slices provide local regular- 
ity conditions for LLG in terms of certain scaling invariant space-time norms of 
Morrey-type, see [551 [Ml • This auxiliary approach leads to partial regularity 
results parallel to those for the harmonic map heat flow in [3T1 [7] , but it fails to 
apply in higher dimensions n > 4. Hence the investigation of more flexible and con- 
ceptually more adapted approaches would be desirable. In this note we establish 
a relationship between LLG and covariant complex Ginzburg-Landau equations. 
This observation is inspired by recent developments in the context of Schrodinger 
maps. It enables us to prove global existence, uniqueness and regularity in arbitrary 
dimensions n > 3 under a smallness condition of initial data with respect to the 
scaling invariant homogoneneous VF^'"-Sobolev norm. Here is our main results: 

Theorem. Suppose A > and n > 3. Then there exist constants p > and c > 
with the following property: Given rrioo G §^ and initial data mo : R" §^ such 
that 

mo~mooe H^n VF1'"(R";R3) 

and such that 

llVmolli-^ < p, 
then there exists a global smooth solution 

m : (0,cx.) X M" 

for the Landau- Lifshitz- Gilbert equation ([1} with the properties that 

sup \\m{t) - moo||//i < ||mo - moo||iji 
t>o 

and 

sup\/t ||Vm(i)||L=o +sup||Vm(t)||L. < c ||Vmo||L" 

i>0 t>0 

and such that 

lim (m(i) - mo) = strongly m n W^i'"(R"; R^). 
The solution is unique in its class and satifies lim m(t) = rrioc in C"'^(R"; R'^). 

// in addition, mo — nioo G iJ°'(M";R'^) for some integer a > ^ + 1, then 

m-moo e C° ([0,oo);i7'"(R";M^)) nC° ((0, oo); i/°°(R"; R^)) . 

Concerning local solutions, existence and uniqueness for if^-regular initial data 
(with a as in the Theorem) has been shown in [SI [T2] in more generality. Observe 
that in particular _ff'^(R") ^ Lip(R"), and local Lipschitz bounds play indeed 
a crucial role in constructing classical solutions for strongly parabolic systems of 
that type. In the context of heat flows local Lipschitz bounds are provided by the 
Bochner identity in conjunction with the maximum principle, and it is possible to 
prove local existence and uniqueness results under weaker assumptions. 
The main point in extending smooth solutions or proving global regularity of weak 
solutions is to show that the smallness condition in VF^'"(R"-) implies a Lipschitz 
bound. In the context of heat flows into an arbitrary compact Riemannian manifold. 
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a global regularity result similar to our Theorem has been proven in |31| Theorem 
7.1 assuming a global Lipschitz bound and a corresponding smallness in energy 
condition on initial data. 

Recall that smoothness and a separate smallness in energy condition on initial data 
alone is not enough to ensure global regularity, as has been shown for LLG at least 
in dimension n = 3, 4 in [S]. In the context of the Theorem it would be interesting 
to investigate whether a bound on the gradient Vm in i°°((0, T); L"(M")) implies 
regularity of an arbitrary weak solution m. In case of heat flows this is a recent 
result, see [36], and motivated by parallel results for the Navier-Stokes equation, 
see [TU]. The argument, however, crucially rests on specific features of the heat flow 
equation. 

In the case of Schrodinger maps, i.e. for A = 0, the parabolicity along with its 
characterisitc local smoothing properties degenerates. Recent global existence and 
regularity results that rest on extensive arguments from Fourier analysis hold true 
under a smallness condition in terms of critical Sobolev norm ij^(R"), see [2113], 
rather than the weaker W^'"{M.^) norm. Thus, in contrast to the present result for 
LLG, results for Schrodinger maps in higher dimensions take into account higher 
order (fractional) derivatives, similar to the case of wave maps, see [^15^ . 
The crucial point and common feature in the proof for LLG and for Schrodinger 
maps is a canonical choice of coordinates on the tangent bundle (rather than on 
the target manifold itself) called the method of moving frames. This method has 
several advantages. First, it reveals the essentially cubic structure of the geomet- 
ric nonlinearity. More importantly in the context of LLG, however, it provides a 
linearization of the quasilinear hybrid-structure and turns the equation into a semi- 
linear (though nonlocal) complex Ginzburg-Landau equation. Finally, the method 
is global, and no small deviation assumption as for the stereographic model is re- 
quired. The strategy of proof is to start from smooth initial data and local classical 
solutions that can be constructed by means of general methods for quasilinear par- 
abolic systems. Crucial higher order energy estimates will be reviewed in Section 
1, where we also address the issue of uniqueness. In Section 2 we shall briefly dis- 
cuss the concept of moving frames. We derive a covariant version of LLG giving 
rise to a covariant complex Ginzburg-Landau equation in terms of a complex vec- 
tor field u that represents the spatial gradient of m and a space-time connection 
field a = (ao,a) that represents a covariant derivative. Using the Coulomb gauge 
div a = we shall deduce some preliminary estimates for the nonlinearities coming 
out of this process. In Section 3 we shall derive corresponding a priori estimate for 
u{t) in the L"(]R") and L°°(R") norm that hold true under a smallness condition 
on ||u(0)||lii. Whereas a priori estimates for Schodinger maps mainly rely on so- 
phisticated Paley-Littlewood decompositions and Strichartz estimates for dispersive 
equation, see [21 3], a priori estimates for LLG will rest on methods from semilinear 
parabolic systems and in particular the usage of weighted- in-time Lebesgue-Sobolev 
spaces as introduced by Fujita and Kato in the context of the Navier-Stokes equa- 
tion, see [121 [H] ■ Accordingly, we obtain in the first instance uniform estimates for 
quantities t >->■ t^\\u{t)\\^2^ for certain S < 1 and 1 1-^ ^/t\\Vu{t)\\Ln featuring the 
appropriate scaling in space-time but improved regularity in space. By Duhamel's 
principle, these quantities can be controlled by ||it(0)||Ln, and we use them to obtain 

Vt\\u{t)\\L^ + \\u{t)\\L'^ < c||u(0)||l,. 

with a universal constant c and independently of t. Transferring this estimate 
back into LLG in Section 4 provides in particular a Lipschitz bound for m, which 
enables us to extend the local classical solutions m uniquely to all times and, by 
approximation, to prove regularity of weak solutions. 
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1. Local existence and uniqueness of classical solutions 

Local smooth solutions for LLG have been used in [5] with the aim of construct- 
ing Schrodinger maps in the limit A \ 0. More recently, an alternative approach 
based on fourth order parabolic approximation has been presented in |19| . The 
approaches work in the context of more general target manifolds equipped with a 
Kahler structure. Here we shall restrict to the case of §^ valued maps where local 
classical solutions can be constructed extrinsically by means of standard methods 
for quasilinear parabolic systems. Our main point is the fact that classical solutions 
persist as long as spatial gradients remain bounded. 

We shall fix constant values for m near infinity, i.e. m(oo) = rrioo for some fixed 
■HT'oo G For this purpose we define, for ct e N, the complete metric spaces 

where _ff'^(R") = Ty^^^(]R") the usual Sobolev space. The distance between two 
elemensts mi and m-z in i/^(M") is given by ||mi — m2\\H^- We also define the 
space 

which is modeled on the Frechet space H°°{W"-) equipped with the usual metric. 



For initial data mo G and a sufficiently large we shall briefly discuss short time 
and global existence and uniqueness of classical solutions. In view of the vector 
identities 

m X Am ~ V • (m x Vm) 

and 

— m X m X Am = Am + |Vm|^m 
valid for m e C^(M";S^), equation ([T]) can cquivalcntly be written as 

(2) — = A (Am + |Vmpm) - m x Am 

which is a quasilinear parabolic system. Taking into account 

m X Am = V • (m x Vm) 

we observe that its principle part is in divergence form and induces a concept of 
weak solutions: 

Definition. We saym e L°°{{0,T); HliR";^'^)) with ^ £ L^{{0,T); L^{W';m^)) 
is a weak solutions of Landau- Lifshitz- Gilbert in (0,T) if 



I dm 

\dt' 



$\ + A(Vm, V$)i2 = (m X Vm, V<I>)l2 + A(|Vmpm,<I>)j 



L2 



for every $ G i/^ n L°°(R";R3) and almost every t G (0,r). 

It is useful to take into account several variants of LLG. Vector multiplication of 
(2]) by m from the left and adding a multiple of ^ yields 

/„% , dm dm , ,9^ ,„ ,9 N 

(3) A — + m X — = (1 + A^) (Am+ |Vmpm) 

from which one can easily deduce the energy (conservation) law 



(4) (i + A2)li?(m(t)) + A£ 



dm{t) 



dt 



2 



= 
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valid for sufficiently regular solutions, hence 

r* r Bm. 2 

dx ds 



(5) E{m{t)) - E{m{0)) ^ 



dt 



1 + A2 

The corresponding a priori estimates are sufficient for constructing global weak solu- 
tions, for which ([5]) turns into an inequality, see [T]. Higher order energy estimates, 
however, are necessary for proving short time existence of classical solutions: 

Proposition 1. Suppose a > ^ + 1 is an integer and niQ G -ff^(K";§2). Then 
there exists a time T > and a classical solution m : (0, T) x M" — > of the 
Landau- Lifshitz- Gilbert equation ([I} with initial data mo such that 

meC"([0,T);iJ,-(R")) and ^ e C° ([0, T); iJ--2(M")) . 
The solution persists as long as ||Vm(t)||Loo remains finite. 

Remark. The terminal time T can be bounded below in terms of WrriQ — mooUff^o , 
where ctq is the smallest integer greater than § + 1- 

The result holds true for a — oo. For less regular initial maps mo e with 
a > ^ + 1 we still have solutions 

me C°((0,T);Fr(R")) 
in particluar, solution from Proposition]^ are in C°°((0,r) x M"). 

A general strategy for solving quasilinear systems by means of spectral truncation 
and higher order energy estimates has been presented in [33| §7 or [S^ Chapter 15 
§7, respectively. Since ([1]) respects the constraint |m| — 1, it will be sufficient to 
construct a solution in a usual Hilbert space by virtue of the ansatz m = §3 + f . 
Obtaining a solution v for some v G C°([0, T); we have to show m — §3 + 1; 

is valued. For this purpose we consider the function 4>{m) = (1 — |mp)2 and 
observe (/'(m) < 2|t;p + jt'l^ which is integrable in space for all t G (0, T). A straight 
forward calculation and integration by parts show 

(j){m{t))dx < 4\\Vm{t)\\loa [ (j){m{t)) dx 

for all t G (0, T), and we deduce from Gronwall's lemma that m is indeed valued. 

1.1. Higher order energy estimates. The proof of persistence given an L°° 
bound on the gradient, is based on the following a priori estimate which can be 
seen as a higher order version of @ . 

Lemma 1. Given any integers a > 2 there exists a constant c > that only depends 
on A and a with the following property: If m G ((0, T); _ff^(M")) is a solution 
of ([D), then 

(6) ^ ||Vm(t)||^„_i + A \\Wm{t)\\l^ < c (l + \\Wm{t)\\l^) || Vm(t)||^„_, 

hence with C{t) — c (l + || VTn(s)|||ac) ds we have for all t e (0,r) 
Jo 

(7) ||Vm(t)||^„_, <e^(*) ||Vm(0)||^„_.. 
Proof. Given any multi-index l<|Q;|<cr— Iwe have 

d°'{m X Vm) ^ m X d^Vm + R where ||i?||i2 < c || Vm||L°° || Vm||^,j-i . 

In fact, the commutator term R only contains mixed derivatives and no absolute 
term m, and we can use the following interpolation inequality (see [34], p. 9) 

\\d^fd-'gh^<c\\f\\^^\\g\\He + \\gh^\\f\\H^ for \/3\ + \j\ = £. 



6 



C. MELCHER 



Moreover, we obtain by a similar argument with a generic constant c 
||<9" (IVmpm) IU2 < c||m||L-||Vm||ioo||Vm|lH-i + ||Vm||i^||Vm||a. 



< c (1 + ||Vm||ioo) ||Vm 



Applying 9" to ([IJ and multiplying it by d°'m we obtain upon integration by parts 

^||9"m||2, +2A||a"Vm||^2 < c || Vm||L- 1| Vm||a.-i ||9"Vm||^2 

+c(l + ||Vm||i..)||Vm||^._i. 

Then ([B]) follows from Young's inequality and after summing over all 1 < |a| < cr— 1, 
while ([7]) follows subsequently from Gronwall's inequality. □ 

A modified version of the Lemma holds true for approximate solutions obtained in 
a suitable contruction process. For a as in Proposition [T] the estimate provides, by 
vfrtue of Sobolev embedding H'^~^{W^) ^ i°°(]R"), the requisite a priori estimate 
for proving local existence. The bound on ||Vm(t)||^„ dt obtained for © upon 
integration in time and the uniqueness result below can be used to bootstrap higher 
regularity estimates in the interior. We refer to [33l [34] for details. 

1.2. Uniqueness and stability. 

Lemma 2. There exists a constant c with the following property: If mi and m2 
are classical solutions as in Proposition^^ on a common time interval [0,r), then 

\\{mi - m2)(i)||^2 < exp c V / || Vm,(s)|j|=e ds \\{mi - m2)(0)||^2 

\ ^=l,2•^0 j 

for every t € (0,r). In particular, solutions in Proposition]^ are unique. If n > 3 
there exists a constant 77 > with the following property: If mi and m2 are weak 
solutions in the sense our Definition in the space L°°((0,T);W^'"{M")) with 

sup |jVmi(i)||L" < 77 
fe(o,T) 



for i = 1,2, then 

\\{mi - m2)(t)||^2 + 2 l|V(mi - m2)(s)||^2 ds < \\{mi - m2)(0)||^2 

for every t e {0,T). In particular, weak solutions in this class are unique. 
Proof. Observe that by ^ and with the notation $ = mi — m2 we have 
$t = A (A<I> + |Vmip$ + m2V$ : V(mi + m2)) - V • ($ x Vmi + m2 x V$). 
Using 2$ as a test function yields, for every t £ (0,r) 

||$(i)||i2 +2A /" ||V$||i2ds < ||$(0)||i2 +2A /" ||$Vmi||2,ds 
Jo Jo 

+2 (^A||$V(mi + m2)||L2 + 11$ x Vmi||L2) || V$||l2 ds 

Now the first claim follows easily from Young's and Gronwall's inequality. 

For the second claim we also use the Sobolev inequality ||$|| < c||V<l'||i2 in 

order to replace the L°° bounds by L" bounds on Vm^ leading to 

ii^wiii^ + A r II v<i>iii2 ds < mmi. + + T?^) c f wvni. ds 

Jo Jo 
for some constant c that only depends on A and n. Choosing rj sufficiently small, 
the second term on the right can be adsorbed, and the claim follows. □ 
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The Lemma is useful in approximating local solutions, and for the time being 
we shall assume mo G i/^(]R";§^). Our aim is to show that classical solutions 
m G C°{[0,T);H-^{R")) are global provided Vmo is sufficiently small in the L" 
norm. A well-konwn approximation device due to Schoen and Uhlenbeck enables 
us to extend this result to weak initial data. 

2. Derivation of the covariant Landau-Lifshitz-Gilbert system 

The formulation of harmonic flows using orthonormal frames on the tangent bun- 
dle of the target manifold under a suitable gauge has been successfully used in 
the context of wave maps (see e.g. [TTl [21]) and more recently in the context of 
Schrodinger maps (see e.g. [371 EH El SI [2S] ) ■ The construction of what we shall call 
the covariant Landau-Lifshitz-Gilbert system is extrinsic and analog to the case of 
Schrodinger maps into . We shall keep our presentation brief and particularly re- 
fer to the self-contained presentation in [3^, Section 2. We shall see that in these new 
coordinates the derivative function will solve a complex Ginzburg-Landau equation. 

Notation. Latin indices k,£ will be used for spatial components. V 

and A will denote the spatial gradient and Laplace operator, respectively. For a G 
R" and u G C" we set a ■ u = '^kUk G C, and accordingly divu = V • u will 

denote the euclidean divergence. Greek indices a, 13 € {0, . . . ,n} will be used for 
space-time components, where a ^ is the time index. Complex space-time vectors 
are denoted by bold letters u = {uq,u) G C"'^^". Finally we set da = gf-- 

2.1. Moving frames. We consider 

mG C°([0,T];i?,°°(R")) and ^ G C° ([0, T]; i?°°(K")) . 

By virtue of a topological construction and a regularization argument, it has been 
shown in f3] §2 that there exist smooth tangent vector fields along m 

X,Y : (0,T) X R" -> T^S^ 

such that 

\X\ = |y| = 1 and X xY, 

that is, the pair {X, Y} forms an orthonormal tangent frame along m. Moreover, 
the vector fields X and Y inherit the regularity properties of m, i.e. 

aaX,a„y G C°([0,r];H°°(M")) for ae{0,...,n} 

We define the map 

a= (ao,a) G C° ([0, T]; i/°°(M"; M1+")) 
with coefficients defined by 

aa^{daX,Y)^~{daY,X), 
where (•, •) is the euclidean product on R"^, giving rise to a unitary connection 

F>a — da-\- iOa. 

Representing the space-time gradient of m in terms oi X -{-iY we obtain coefficients 
(8) Ua — {dam,X) + i {dam,Y) 

such that 

dam = Re{ua) X + lm{ua) Y 
The coefficients Ua form a smooth complex vector function 

u= (uo,m) G C° ([0,T];i7°°(M";Ci+")) . 
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Taking into account the relations 

daX = —Ke{ua) m + QaY and daY = —lm{ua) m — a^X 
we obtain the zero torsion identity 

(9) DaUp = DpUa 

and the curvature identity 

(10) Rap := [Da, Dp] = i [daap - dpaa) = ilmiuaUfi). 

2.2. The covariant LLG system. Suppose T £ (0,r*) and 

me c"([o,r,);irr(R";S')) 

is a solution of LLG with m(0) = m-o G H^{W;S'^) as in Proposition [H On [0,r] 
the construction of a moving frame {X, Y} along m from Section [^TTl can be carried 
out. A straight forward calculation shows 

n n 

Am = ^ {dkRe{uk) - at lm{uk)) X + {dklm{uk) + ak Re(ufc)) Y + |apm. 
fc=i fc=i 
Upon multiplication hy X + iY , the Landau-Lifshitz-Gilbert equation 
dim 

—— = A (Am + iVmpm) — m x Am 

at ^ ' 

turns into the following system that we shall call the convariant LLG system. 

Proposition 2. Suppose m is a solution of the Landau-Lifshitz-Gilbert equation 
as in Proposition[l\ on an interval (0,T,) with initial data mp G -ff^(M";§^). For 
T € (0,T*) and with the notation from Section \2.1\ 

(u,a) e C° ([0,T];i?°°(R";Ci+" x R^+")) 

solves the covariant Landau-Lifshitz-Gilbert system: 

Mo = (A - i) J2k=i ^kUk 

(11) < DpUa = DaUp 

daop - dpaa = lm{uaup) 
Hence u = (ui, . . . ,m„) solves the covariant complex Ginzhurg- Landau equation 

n 

(12) DQUi^{\-i)Y,{DkDkUi + RaUk) for le{l,...,n} 

k=l 

and attains initial values u{0) — (Vmo,X) +i(Vmo,y) in _ff°°(R";C"). 

2.3. The Coulomb gauge. We observe that the covariant LLG system in invari- 
ant with respect to the individual choice of the moving orthonormal frame {X, Y}. 
This local rotation invariance translates into the fact that (jlip is invariant with 
respect to gauge transformations 

u I— ?> e~*^M and i— > + dad. 

A canonical choice is the Coulomb gauge which is characterized by the equation 

n 

div a* ~ diUi — 
e=i 

and obtained from our initial moving frame {X, Y} by a gauge change represented 
by a suitable gradient 

V6I G C"([0,T];i?°°(R")) 
where 6 solves, for all t £ [0,r], an elliptic equation 

-Ae{t) ^ div a{t) in R". 
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Note that V9{t) can be obtained by taking double Riesz transforms of a which are 
bounded operators on every Sobolev space i/'^(M";R"), see e.g. OH^. It follows 
that the new frame {X* ,Y*} inherits the regularity of {X, y}, same for a* and u*, 
see [3] Proposition 2.3. In the sequel we will skip the superscript *. 

2.4. LP estimates for a(t). RecaU from dTU]) that for a, /3 e {0, . . . , n} 

daap - dpaa = Im(MQt2^) 
which implies that under the Coulomb gauge for a — . . . ,n 

(13) - Aoq = dklrtv{uaUk) 

A priori estimates in Section [3] take into account higher power norms of u{t), 
i.e. p — n/5 larger than the critical power n, and L" norms of Vu{t). In turn we 
will need corresponding estimates for a{t) in terms of u{t). 

Lemma 3. Suppose 6 S 1) and ti > 3. Then we have under the Coulomb gauge 
with a constant c that only depends on S, A and n 

(14) \\ait)\\^^<c\\uit)\\K 

for all t £ [0,T]. Moreover, there exists a decomposition ag = Oq"'^' + a^'^ where 
e C"([0,r];LT(M")) and e C°([0, T]; i^fe (R»)) 

and such that 

(15) ||aW(t)||^. <ch(t)||^.||V^t)|U. 
and 

(16) \\a^^\t)\\^^ <c\\uit)\\K 
for all t £ [0,T]. 

Proof. All estimates are of elliptic type, and we shall suppress the time dependence. 
Recah that if 1 < p < n and / e LP(R";R") then 

-Au^div/ in M" 

has a (weak) solutions v £ VF^'P(M") given by 



(-A)-^5;]7^fe/fc 



fc=l 



where the TZk are the Riesz transforms and (—A) 2 is the Riesz potential cor- 
responding to the Fourier multiplier 1/|^|. It follows from the Hardy-Littlewood- 
Sobolev inequality, see e.g. [30] p. 119, that v £ L"-p (M") with 

(17) <c||/|Up. 

Moreover, the solution v is unique in its class. Hence, in view of (1131) . 

n 

— Au£ = 9fcIni(i2fcUf) — divlm(uu^) 

and estimate ([T4)) follows from ([T7| with p — For the decomposition of oq we 
observe that by Proposition [2] 

Im(uMo) = lm^{X ~ i) u DiUg^ 

= Im^(A — i) u divu^ + Re^(A — i) {a ■ u) ilj . 
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Accordingly we define and Aq by the equations 

—Aa^Q ^ — div ^AIni(u divu) — Re(M divu) 

and 

—Aa^Q^ = div ^ARe(M (a • u)) + lm(u [a ■ u)) 

Estimate (flSl) follows again from Holder and p7|) where p = j^. Taking into 
account estimate ([TB| follows with p = jgzj- The proof is complete. □ 



Remark. For Schrodinger maps, i.e. for A = 0, we would obtain in view of (|9]) 

—Aoq = div div id — Re(?2 (g) u) 

and hence \\ao\\ LP < c \\u\\ for all p ^ (l,oo) by the C alder on- Zygmund inequality. 
The fact in case of LLG the estimate for oq takes into account Vu is not essential 
for our regularity argument. 

3. Estimates for the covariant Landau-Lifshitz-Gilbert system 
In this section we consider 

M e C° ([0,T];i?°°(]R";C")) 
solving the covariant Ginzburg-Landau system (jl2p under the Coulomb gauge. 

3.1. Linear estimates for the semigroup. Let us first recall the fundamental 
estimates for the dissipative Schrodinger semigroup S — S{t) which is generated by 
(A — i)A. With a slight abuse of notation, we represent this semigroup in terms of 
the Fourier multiplier 

which is a Schwartz function for every A > and t > 0. Thus the associated kernel 
St is also a Schwartz function with the scaling property 

Stix)=t-''/^S{x/Vi). 

Writing S{t)f — St * f for f & iP(M";C) we obtain from Young's inequality the 
following mapping properties on Lebesgue spaces parallel to those of the heat kernel 
(see also [37] Proposition 2.3): 

Lemma 4. Suppose A>0, l<p<(7<cx) and a £ N. Then there exists a 
constant c > such that 

\\'^''S{t)\\c(LP;L'>) <ct-^(p-^y^ for all t > 0. 

3.2. Nonlinear estimates. We recall that the covariant LLG system in Propo- 
sition [2] gives rise to a complex Ginzburg-Landau system (|T2|) which, under the 
Coulomb gauge, can be written as 

(18) ^ = {X -t) An + F{a,u) 

at 

where, for a = (ao,a), the nonhnearity is given by 

(19) Fi{a.,u) (A - i) |i ^ (lm{ueuk) Ukj +2i{a- V)ui - |apuf | - ia^ut. 

Taking into account the decomposition of ao — a^^^ + Oq^^ from Lemma [3] we find 
that F splits into five terms 

F{&,u) =f(i) + --- + F(5). 
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Apart from the cubic term F*^^-' that only depends on u in a local fashion, the 
functions F^^^ to F*^'^-' can be considered as nonlocal multilinear operators acting 
on u, that is a = a(u) is considered as a function of u. More precisely, in view of 
Lemma [21 F^^^ and F*^^^ are quadratic in u and linear in Vu, while F(3) and F^'^) 
are quintic in u. Therefore we introduce the functions 

n 

fc=i 

/(2) ^(2) ^_ ^(4) ^ {X-i)2i{a-V)u-ia'i^^ u, 
/(3) ^(3) _^ j^(5) ^ -(\-t)\a\\~ia^^K 
and obtain from Lemma |3] and Holder's inequality: 

Lemma 5. Suppose n > 3 and S G (i, 1). Then there exists a constant c such that 

(20) \\f'Ht)h^ < c\\uit)\\K, 

(21) \\f^'Ht)h^ < c\\u{t)\\K\\Wu{t)\\L^, 

(22) ll/'''(i)IL^ < c\\u{t)\\K 



for allte [0,T]. 



3.3. Duhamel's principle and Pujita-Kato type estimate. At this point we 
have essentially reduced the problem to a linear one. Writing / = + + 
we obtain the following representation formula 

(23) u{t)^ S{t)u{Q) + S*f{t) 
where S(t)f — St*f (i.e. convolution in space) and 

(24) (^ */)(*):- f S{t^s)f{s)ds 

Jo 

i.e. convolution in space and time. Based on this representation we shall derive a 
priori estimates in suitable (scaling invariant) weighted-in-time Lebesgue-Sobolev 
spaces. For this purpose we set, for (5 e (5, 1) and i S [0,r], 

K{t) = sup t-^\\u{t)\\^^ 
K'{t)^ sup T^|lVu(r)||i. 

re(0,t) 

and 

R{t) = max{K{t),K'{t)}. 
It follows from Sobolev embedding that 

(25) Rit) <ct^\\uit)\\^^ for ah t € [0,T]. 

Hence R : (0,T) [0,oo) is continuous with limi?(i) ~ 0. We also introduce 

Ro{t)^max\ sup r^||5(r)u(0)||^n , sup T^|jV5(T)?/(0)||L. I ■ 

[re(0,t) re(0,t) J 

According to Lemma S] we have 

(26) Roit) <c\\u{0)\\l^ 

independently of i > and for a constant c that only depends on 5, A and n. 



12 



C. MELCHER 



Lemma 6. Suppose S G (|, |). Then there exists a positive constant cq such that 

(27) R{t) < Roit) + Co {R{tf + R{tf) . 

Moreover, there exits a positive constant with the following property: If 

sup Ro{t) < ro 

tG{0,T) 

then R{t) < 2 Ro{t) for every t e [0,T]. 

Proof. Estimate (P7)) will follow from ([25)1 . Taking into account that for a,[3 < 1 
(restricting the range of admissible S) and for positive t with a constant c = c(a, /3) 

{t~s)-"s-^Us = ct^~°'-'^ 

we obtain form Lemma HI and Lemma [S] with a generic constant c 

(5 * (t)||^„ < c if (t)3 - = cif(0' t^, 



(s'*/(2')(t) ^„ <cK{tfK'{t) / (t- = cif(<)2if'(i) 



* /^')) {t) „ < c if(t)5 / - = cif 

V I Jq 

(0,T), provided (5 € (|, 1) (used in the 3rd estimate). Moreover 

V (5 (t)||^^ <cK{tf j (t-s)"¥s-i(i-«) ds = cif(t)3i-^ 

(5 * (t)||^^^ < c K{tfK'{t) f\t - s)-*s-i+^ ds = cK{tfK'{t) t'^ , 



(s * (t) ^ < c i4:(t)^ /" (i - s)-^(^''-2)s-i(i-<5) ^ cK{tf t-5 



for all t e (0,r) provided (5 e (|, |) (used in the 1st estimate). Hence (P7)) follows. 

To prove the second claim, we argue by contradiction and assume there exists 
to e (0,r) such that i?(to) = 2i?o(to) 7^ 0. Then, by virtue of (EH), 

2i?o(t) < i?o(to) + co (8i?o(to)' + 32i?o(to)*)i?o(to)- 

Thus for < i?o(to) < ''o with ro sufficiently small 

1/co < 8 i?o(<o)' + 32 Roitat < 1/co, 

a contradiction. Since i? is continuous with limt\,o Rit) = 0, the claim follows. □ 

Recall that by (^5)) the smallness condition can be expressed in terms of ||u(0)||l»i. 
Hence we obtain: 

Proposition 3. Suppose S G (|,|). Then there exist positive constants p and c 
depending only on 5, A and n with the following property: If 

ll«(0)|k" <P 

then the following estimate holds true for every t e (0, T] 

t^\\u{t)\\^^ +ti\\Vu{t)\\L^ <c\\u{0)\\l^. 
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3.4. Uniform bounds in and Holder continuity. 

Lemma 7. Suppose a G (0; Then there exits a constant p and c depending only 
on a, A and n with the following property: If 

h(0)||L" <P 

then the following estimate holds every t € {0,T] 

t'^ws/um^^ <c||u(o)iu.. 

Proof. We first show that for some constant c > independent of i e (0, T] 

(28) t^||Vu(t)||^^ < c||u(0)||l. +c {R{tf + R{tf) . 

To this end we repeat the argument in Lemma [5] with a change in the exponent of 
t — s according to Lemma S) With a generic constant c we obtain for all t g (0, T] 



7 * /(I)) {t)\\^^ < c R{tf J {t- 5)-^(3*+")s-i(i-*') ds = cR{tf r"-^, 
V [S * {t)\\^_^ < c R{tf J\t - 5)-i(25+«)5-i+5 ds = cR{tf t-^, 
V [S * it)\\^^ < c R{tf At- ds = cRitf r"-^ 



The first estimate requires 3(5 + a < 2 while the third estimate requires i5 > | which 
is possible by an appropriate choice of (5 € (f 1 1) provided a e (0, -1). Finally we 
observe that according to Lemma 0] 

t"-^\\VS{t)um\^^ <c||u(0)|U. 

for every t > Q and with a constant c that only depeneds on a and n. Hence 
Duhamel's formula implies (|28l) . If p is sufficiently small Lemma [S] implies 

R{t)<2Roit) while by (IMl) i?o(i) < c ||w(0)||l" < cp, 

and we obtain II VM(t)||^^ <c {l + p'^ + p^) \\u{0)\\ l'^ as claimed. □ 

The Lemma implies, by virtue of Morrey's inequality, Holder continuity of u(t) for 
every t £ (0,r) with Holder exponent a € (0, ^) and bounds that only depend on 
||u(0)||l'i, a, A, n, and t. Along the lines of the above argument we also obtain 

(29) ||w(t)||L" <c||u(0)||tn for all te[0,T] 

provided ||u(0)||l'> < P- An L°° bound for u{t) can be obtained by scaling. In fact, 

\Ht)\\^^ < \\uit)\\c^ < c (||u(t)||i„ + ||V«(i)||^^) 

hence, see m P-9, h(i)ILoo < e \\u{t)\\J^ \\S/u{t)\\^ < ^||w(0)||i.. 

Proposition 4. There exist constants p and c depending only on A and n with the 
following property: If 

lk(0)||L" <P 

then the following estimate holds true for every t G [0,T] 

Vt ||m(<)||l- + ||w(OIU" < c||u(0)||l.. 
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4. Proof of the Theorem 

4.1. Initial data in iJ^(M";S^). Given a map mo in this class, there exist, ac- 
cording to Proposition [TJ a terminal time T, > and a smooth solution 

meC"([o,r,);i?r(R")) 

with m(0) = mo. For T E {0,T^,) we obtain, by Proposition [21 a solution 

ue C°([0,r];iJ°°(R";C")) 

of the corresponding covariant complex Ginzburg-Landau equation under the Coulomb 
gauge. Since |Vm| = |u| we have for t G [0,T] 

(30) ||Vm(0||Loo < ||u(i)||L~ while \\u{0)\\l'^ < ||Vmo||L- 
Assuming ||Vmo||L" < P, Proposition |4] applies, and we obtain 

(31) ||Vm(i)||^<^ < -^llVmollin for aU te(0,r) 

where c only depends on A and n. Since T G (0,T,) is arbitrary we deduce from 
Proposition [T] that m can be extended to all times. In turn, ([5T|) holds true for any 
t > 0. In conjunction with the energy inequality ([5]) we obtain 

(32) ||m(t) - moo||L~ < c||Vm(t)||^i ||Vm(t)||£2 < ct ||Vmo||L>>nL2 
for all t > 0. Hence m{t) moo in C^(R") as t cx). 

4.2. Initial data in _ffJ(R";§^) for a > ^ + 1. Given a map mo in this class we 
find, according to Proposition [TJ a terminal time > and a local solution 

meC"([0,r,);i/,-(R")). 

In order to obtain a uniform bound for Vm, we shall use the result in Section WTU 
in conjunction with the following approximation result, originally due to Schoen 
and Uhlenbeck, see [28] and [31] Proposition 7.2. 

Lemma 8. Suppose m : R" -)> is suc/i f/iaf m - m^o € i?^ n W^i^"(M"). Then 
there exists a sequence of maps m}'^^ £ i?^(R";§^) such that 

If m £ i/^(M") for some a > I then {m^''^} is uniformly bounded in iJ^(R"). 

Proof. With a standard moUifier ip G C^{Bi{0)) and <y3e(x) — e^^'^if{x/e), we 
obtain, taking into account m — m^o G L^(M"), 

* m — moo = * {m — moo) G iJ°°(M"). 

Repeating the argument from [33] yields for all x £ R" 

(33) (l-|(^e*rn)(x)|)"=dist((^, *m)(a;),§2)"<c /" |Vm|"dx^O 

uniformly as e — > 0. For e sufficiently small we have inf jiyS^ * m| > i and let 

mW:=^l^ for e = \ 0. 

We *m| 

We observe that we can represent m*^''^ — moo = ^ifsk * — moo) with a smooth 
function F : R^ — > R^ from which we can assume -F(O) = and which induces a 
bounded mapping on L°°n-ff'^(R"; R'^) for every cr e Z, see e.g. [34] p. 11. Moreover, 

Im^*^-* — m| < jm'-'^'-' — ipe * m| + * m — m| < 2 * m — m|. 
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hence m^^^ —m — > in L^(R") as A; — > oo. Finally, |V [ipe * 'm) \ > |(pe*m| \\'m^'''>\ 
and therefore by ([55]) we obtain, for 2 < p < n, 

limsup llVm'-'^-'llLP < lim ||iy9e^ * m||ip = ||m||ip < liminf || Vm^'^-' ||iP. 

^—^OO /i'— 7-00 fc— >oo 

It follows that m^'^) - m ^ in iJ^ n VFi'"(R") as fc ^ oo as claimed. □ 

If we assume || Vmo||^„ < p with p as in Proposition^ then, for large k, the same is 
true for a approximating sequence mg'^'' of initial maps as in Lemma [S) We obtain 
a corresponding sequence of global solutions m'--''^ e oo); i/^(R")) such that 

||VmW(i)||Loc < -2=||Vm{,'')||L" for ah t > 

and a universal constant c. Passing to the limit fc — > oo we see that by Lemma [2] 

m('=)^m in L^^{{0,n); L^R"")) 
and by lower semicontinuity of norms 

sup ||Vm(t)|l^^ < oo. 

We infer that m persists for all times by virtue of Proposition [1] and is unique in 
its class according to Lemma Similarly, (PT|) and p2p carry over to the limit. 

4.3. Weak initial data. Given intitial data mo : R" such that 

mo-rriooeH^n iyi'"(M";R3) and HVmoU^,. < p 

with p as in Proposition 21 there exists, by virtue of Lemma [51 a sequence of 
approximating initial maps rrifl'^ G H^{M") and, for large fc, corresponding global 
solutions € C°{[0,oo); H^{W)). If cp < 1] with 7] as in Lemma [5] we deduce 

in conjunction with the energy inequality ([5]) 

\\m^''\t) - moollffi < Wm'^o^^ - rriooWm 

and 

""^ amW... ' . ^ 1 + , (fe). 

dt < — - — E(m\^ ') 



-it) 



L2 



dt 

while from Proposition [H 

Vi||Vm«(t)|U=o + ||VmW(t)|U,. <c||Vmf |U„ 

and from 

\\m^''\t) - mooIlL- < ct-^ W^m'-Q ^WL^nL^ 

for all i > and with a univeral constant c. Now it is easy to see that any weak 
limit m is a global weak solution of ([2]) such that the prescribed initial data tuq 
is continuusly attained in and such that all the above estimates are satisfied by 
m. Hence, by Lemma [2l the solutions m is unique in its class. 
We aim to prove uniform local bounds in spaces of space-time Holder continuous 
functions. For this purpose we invoke the following localized energy inequality for 
Landau-Lifshitz-Gilbert, that can be found in [24; Lemma 2 or [3S] Lemma 5.10: 

Lemma 9. Suppose m is a smooth solution of the Landau-Lifshitz-Gilbert equation 
in a space-time cylinder 

Pr{zo) = [to, to + r^) X Brixo) where zq = [ta.xo) 

then, for a constant c that only depends on X and n, 



(34) / |Vm|^ dx 



dm 



dt 



2 



C 



dz < — I iVml^ dz 



Prize) 
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Given any compact subset Q C (0, oo) x M" there exists a constant c that only 
depends on || Vm''^) ||ioo(g) such that, for ah Pr{zo) C Q 



and hence, by the parabohc version of Morrey's lemma, see IB], a uniform Holder 
bound of {m'^'^^} locally in Q. Hence m is locally Holder continuous with locally 
bounded gradient. But then the linear theory of parabolic systems as e.g. in [22] 
and a bootstrap argument implies m is smooth in (0, oo) x R". 
Finally, we show that initial data is strongly attained in n Indeed, 
by lower semicontinuity of norms and the strong convergence of initial maps 

||mo — mooll < limsup ||m(t) — mooll < limsupliminf ||m'''^''(t) — mooll 
< limsup ||m[)''-' - mooll = ||mo - mooll 

k—^cx) 

where || ■ || is the H W^'"^^ norm. This completes the proof of the Theorem. 



|VmW|2+r2 



dt 



dz < cr 



2+n 
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